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The propagation of ultra-short optical pulses in nonlinear media is of great importance in 

using femto-second laser pulses in air [3|, U, 15], and inertial confinement fusion [6j. The 
nonlinear propagation of ultra-short optical pulses can be modeled by using the nonlinear 
Schrodinger equation with short-pulse nonlinear derivative corrections, which can give rise to 
the filamentation of light pulses, and subsequent formation of light pipes. Modulational and 
filamentational instabilities therefore sometimes pose problems concerning pulse propagation 
in nonlinear media. Thus, for propagation times much longer than the typical pulse length, 
it is of importance to find means for the nonlinear stabilization of such optical pulses. 

In this Brief Report, we investigate the statistical properties of a generalized nonlinear 
Schrodinger equation, taking into account the Kerr nonlinearity, linear absorption, nonlinear 
dispersion, delay in nonlinear refractive index, third-order dispersion, differential absorption, 
and diffraction. The equation is of relevance for Raman excited four-state atomic systems. 
We derive an equivalent wave-kinetic equation that governs the propagation of optical quasi- 
particles, enabling us to study the effect of partial pulse coherence. The newly derived 
equation is analyzed for the case of spectral pulse broadening, and it is found that the 
relevant growth rate can be significantly reduced if an appropriate random phase techniques 
are used. 

Recently, Hang et al. [3] derived a generalized nonlinear Schrodinger equation (NLSE) 
which describes the nonlinear propagation of optical pulses in Raman excited four-state 
atomic systems (see Fig. 1). In a dimensionless form, the NLSE can be written as 

id z u + d 2 u + 2\u\ 2 u = —i [d u — did t (\u\ 2 u) — d2ud t \u\ 2 — g^c^m] + d±d t u — d^V\u, (1) 

where u is the normalized pulse amplitude and dj for j = to 5 is the dimensionless 
coefficients for linear absorption, nonlinear dispersion, delay in nonlinear refractive index, 
third-order dispersion, differential absorption, and diffraction, respectively. Equation (JT]) 
thus describes the nonlinear propagation of ultra-short optical pulses when the terms in the 
right-hand side of (1) are significant (i.e. dj ~ 1) 16]. In what follows, we will look at the 
one-dimensional problem, and thus set d§ = 0. 

In order to analyze the spectral evolution of partially coherent optical pulses, we next in- 
troduce the Fourier transform of the two-point correlation function, i.e. the Wigner function 
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FIG. 1: Energy-level diagram showing the excitation of the relevant four-state atomic system. Here 
Up is the center frequency of the weak probe field, while ujb and u>c are the frequencies of two strong 
cw control fields. The one-photon detuning frequencies are defined according to Auj p = cu^i — uj p 
and Auj-q = <jj\i — wb- Moreover, 2f2 p , 2S7b, and 2f^c are the Rabi frequencies for the indicated 
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transitions, and two-photon resonance is assumed to always be maintained 
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where the angular bracket denotes the ensemble average M . The Wigner function represents 
a generalized distribution function for optical quasi-particles. The Wigner method |13| . 
as well as the equivalent mutual coherence method [14j ]. have been used to analyze the 
modulational instability in the "standard" nonlinear Schrodinger picture, in which the right- 
hand side of Eq. ([1]) is identically zero. From the definition (j2]) follows the normalization 

I(z,t) = (\u(z,t)\ 2 ) = J dup{z,t,u), (3) 

giving the pulse intensity in terms of the Wigner function. 

Applying the z-derivative to the definition ([2]), we obtain the wave- kinetic equation for 
optical quasi-particles 

d z p - 2u)d t p - 4/ sin (\ dtd<A P = ~2d p - d 1 2ujI sin ( \ dtdJ] P - d t (/cos ( \ dtd<A P 



+2d 2 {d t I) cos (I dtdJ) P~d 3 (3u 2 - \d 2 t ) d t p - 2d 4 up. 



(4) 
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The wave-kinetic equation (j4j) determines the phase space evolution of partially coherent 
optical wavepackets. In the low- frequency limit, we can retain only the first term in the 
operator expansions to obtain 

d z p - 2ud t p - 2d t I d w p = -2d p - d\ [ud t I d^p - d t (Ip)} + 2d 2 pd t I - 3d 3 u 2 d t p - 2d 4 up.(5) 

Here we see that the terms in the left-hand side of (]5]) resemble the terms in a Vlasov 
equation. Moreover, Eqs. (jlj) and (jSJ) have a time-independent solution 

p(z,u) = p {u)e- 2{do+d4Uj)z , (6) 

exhibiting a spatial diffusive influence of the linear and differential absorption. For a coherent 
distribution, i.e. po(u) = Iq8{lo — uj q ) for some frequency ljq, we obtain the intensity I(z) = 
Joexpf— 2(d + diUJo)z\, while for a Gaussian distribution po(u) = (Io/V^rr A) exp[— (cu — 
a>o) 2 /2A 2 ] with the spectral width A, we obtain the intensity I(z) = Iq exp[~2(d + d^uj — 
d±£*?z)z\. We note that the finite spectral width in the latter case gives rise to growing 
behavior for after a certain distance of propagation. In fact, the spatial intensity decay halts 
for z cr it = (do + c?4Co>o) / d^A 2 , indicating the breakdown of the above model, as we do not 
expect growing modes from loss terms. For a broad spectral distribution, this distance may 
be short. 

Next, we analyze the modulational instability of broadband optical pulses that are gov- 
erned by Eq. (]4]). Letting p(z,t,u>) = p(z,u) + pi(z,u) exp(iKz — iQt) in the latter, where 
I Pi I <C p, we linearize the wave kinetic equation against p\. If the dimensionless perturbation 
wavenumber satisfies K 3> do, d&u, the perturbation wavelength is much smaller than the 
decay length. Thus, we may treat the background distribution as a constant, and neglect 
the terms containing do and d 4 u. We then obtain the nonlinear dispersion relation 

~ J ^ K + 2uQ + d^U - d 3 (3a; 2 + VL 2 /A) VL ' 

(7) 

The case of higher order dispersive corrections to the nonlinear Schrodinger equation in the 
Wigner picture has been treated in the wavenumber domain in Ref. 151 ] . 

In the case of coherent optical pulses, we let Po(oj) = Io5(u — <^o)- Here the dispersion 
relation (j7j) reduces to 

K = — [2cu + (2di + d 2 )I - d 3 (3u^ + fi 2 )] SI 

± [(di + d 2 ) 2 I 2 Vt 2 + (1 - 3d 3 ^o) 2 ^ 4 - 2(2 - diw )(l - 3d 3 u )m 2 } 1/2 . (8) 
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For simplicity, we focus on the case u> = 0, at which we have the growth rate T = — Im(i^) 

r = [Am 2 -n 4 - + d 2 ) 2 i 2 n 2 ] 1/2 , (9) 

where the nonlinear dispersion and delay in the nonlinear refractive index gives a stabilizing 
effect to the regular Kerr modulational instability (see Fig. 2). We note that the modula- 
tional instability follows as a consequence of the nonlinear pulse propagation equation (1) 
governed by the underlying dyamics of the four-state atomic system [8[ . The latter system is 
in the form of a four-wave mixing set of equations, and it is therefore perhaps not surprising 
that the modulational instability occurs for the corresponding nonlinear pulse propagation. 
However, other level configurations may in principle posses similar instability properties, but 
the issue of the uniqueness of such instabilities to the nonlinear four-state atomic system is 
left for future research. 

If the background pulse \l/o has a random phase, with a coherence width A, this corre- 
sponds to a Lorentz distribution 

For the case of = 0, i.e. vanishing differential absorption, we can integrate the dispersion 
relation (JTj) for po given by the Lorentz distribution (TTOT) to obtain 

K = - [(2^ + d 2 )I - d 3 tt 2 ] - 3d 3 A 2 tt + 2iAtt 

± [(di + d 2 ) 2 I 2 Q 2 + (1 + 3id 3 A) 2 fl 4 - 2(2 + zd x A)(l + 3id 3 A)IQ 2 ] 1/2 . (11) 



From the dispersion relation ffTTj) we see that the spectral broadening interact in nontrivial 
ways with the modifications of the nonlinear Schodinger equation. For a pure Kerr nonlin- 
earity, i.e. dj = for all j, the random phase of the b ackg round pulse will give rise to a 



reduction of the modulational instability growth rate [13l,ll4j]. Here we find that the random 
phase will also contribute to the real part of K, and the imaginary contribution via interac- 
ions between the corrections to the Kerr nonlinearity and spectral broadening. As noted in 
I5I , higher order dispersive effects may affect the modulational instability growth rate, but 
on the perturbative level, third order dispersion only gives a shift in the real wavenumber. 
This can be seen by setting di = d 2 = 0, and d 3 A, d 3 Q -C 1. Linearizing Eq. (11 II) we obtain 

0,7 An 2 

K = d 3 Q 3 - 3d 3 A 2 Q + ^AQ ± i (AIQ 2 - ^f' 2 ± 2( ^ 2 _^ 4)1/2 (™ 2 - 42) , (12) 
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FIG. 2: The coherent modulational instability growth rate T plotted as a function of Q, using Eq. 
(jlip with A = (see Eq. ([9])). The full thick curve has d\ = d,2 = d% = = 0, giving the Kerr 
modulational instability growth rate, while the dashed curve has d\ = d% = 1 and d 3 = d^ = 0. We 
see that the effect of nonlinear dispersion and delay in nonlinear refractive index is to reduce the 
modulational instability growth rate. 

and we see that the third order dispersion only contributes to the real wavenumber shift. 

We have solved for the modulational instability growth rate numerically and plotted the 
result for different parameter values in Figs. 2 and 3. We note that when A = we regain 
the coherent growth rate (Q. Moreover, the third order dispersion only couples to A and is 
thus not present in the coherent case. 

From Fig. 3 it is clear that for d 3 = 0, the partial coherence of the background pulse 
gives rise to a reduced growth rate, and can thus act as a means for stabilizing optical pulses 
in four-state atomic systems described by Eq. ([1]). The effect of third-order dispersion is 
to further reduce the instability growth rate, as expected, and the new branches in the 
modulational instability growth that occurs for high frequency perturbations, i.e. Q > 1, are 
not valid, as the assumptions underlying the derivation of Eq. ([I]) are no longer satisfied. 
Thus, the strongly growing modes depicted in Fig. 3 for Q > 1 are not physical. 
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FIG. 3: The incoherent modulational instability growth rate T plotted as a function of fi, using Eq. 
(jlip . for different parameter combinations. The full thick curve is the coherent Kerr modulational 
instability growth rate, i.e. d\ = c?2 = cfe = d± = A = 0. The remaining curves has A = 0.1 and 
c?4 = 0. From the top the dashed-dotted curve has d\ = efo = and d% = 1, the dashed curve has 
d\ = c?2 = c?3 = c?4 = 0, the dotted curve has ^1 = ^2 = ^3 = 1, and the curve with long dashes has 
d\ = d2 = 1 and d^ = 0. The strongly growing modes for O > 1 are however not physical, as the 
underlying assumptions of Eq. ([1]) are no longer satisfied [see Eq. (|12p ] . 

In conclusion, we have presented an investigation of the modulational instability of broad- 
band optical pulses in Raman excited four-state atomic systems. For this purpose we have 
obtained a wave-kinetic equation from the modified NLSE, which accounts for numerous 
non-ideal effects embedded in the right-hand side of (1). Using standard technique, we then 
derive a nonlinear dispersion relation from the wave kinetic equation. The nonlinear dis- 
persion is analyzed for coherent and broadband spectra of optical pulses. It is found that 
the growth rate of the modulational instability is reduced in the presence of a Lorentzian 
optical pulse distribution. A reduced growth rate insures stability of optical pulses over long 
distances in four-state atomic systems. 
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